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ABSTRACT 


The number of nonlinear loads connected to power distribution systems is rapidly 
increasing. These loads disturb the sinusoidal nature of the power distribution systems 
to which they are connected. Consequently, all connected loads may be affected by 
even a single nonlinear load. This makes power quality an increasingly important is- 
sue, particularly in shipboard applications because equipment malfunction can cause 
the loss of life. 

This study primarily develops a modeling technique for use with the harmonic 
power flow program (HARMFLO) developed by the Electric Power Research Institute. 
A land based power distribution system is modeled and simulation results are com- 
pared to field measurements. Additionally, the causes and effects of harmonic distur- 
bances in power distribution systems are reviewed. 

In spite of the fact that many parameters of the system being modeled are not pre- 
cisely known, simulation results are similar to the field measurements. This implies 
that HARMFLO has applications in isolating existing harmonic disturbances of installed 


systems and impact assessment for planned systems with one or more nonlinear loads. 
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I. INTRODUCTION 


A. OVERVIEW 

All ships and aircraft in the United States Navy have ac power distribution systems 
installed. Each of these systems is subjected to loads that produce disturbances in the 
system current flow. Some of these disturbances are referred to as harmonic currents 
and can cause equipment malfunction, power loss, and even catastrophic equipment 
failure [1]. Therefore, a tool is needed to assess the impact of harmonic producing 
loads on power distribution systems. 

The Electric Power Research Institute (EPRI) produced a harmonic flow computer 
program called HARMFLO. Power distribution systems can be modeled and simulated 
with this program. Since its development in the early 1980's, HARMFLO has not been 
exhaustively tested and compared to field data on power distribution systems [2:pp. 1— 
2}. Asa result, the effects of program limitations and accuracy are not quantitatively 
known. If the accuracy of HARMFLO is documented, program simulations may be used 
to isolate existing problems with installed equipment and to estimate the impact of 
proposed harmonic producing loads on power distribution systems [3]. 

This study develops a power distribution system modeling technique for use with 
HARMFLO. A land based power distribution system is modeled in detail and exten- 
sively simulated. Many details about the system are not known precisely; much of the 
modeling process involves making assumptions and estimating system parameters. 
This lack of specific system information is consistent with other studies of power dis- 
tribution systems. Simulation results are compared to field measurements and conclu- 


sions drawn about program accuracy and limitations. Additionally, HARMFLO output is 


compared with current military power quality guidance; comments about the applica- 


bility of such comparisons are provided. 


B. THE FIELD STUDY 
1. The System 

The system under test is the Union Electric Audrain County 12 kV distribution 
system. A number of residential, commercial, and industrial loads are connected to the 
distribution system. One industrial load, hereafter referred to as the customer, operates 
seven 100 HP adjustable speed DC drives. Additionally, there are six capacitor banks 
installed throughout the system. A one-line drawing of the Audrain County subsystem 
is shown in Fig. 1. [4] 

The system connections shown in Fig. 1 are modeled as transmission lines. 
Impedances are taken from the utility database and are expressed as a percentage of the 
system impedance base [4]. The system impedance base is defined as 


Z 
zane) (1.1) 


S530) 


where Vay), is the base line to line voltage in kV and S,.,9) is the three phase system 


power base in MVA [5:pp. 32]. The system power base is 100 MVA with a voltage 
base of 12.5 kV [4]. Therefore, the system impedance base is 


Pa 
PRP ha ane (1.2) 
100MVA 


and the system current base is 
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Fig. 1: One-line drawing of the Union Electric Audrain County 12 kV distribution system 


The loads are estimated from utility data and are based on an assumed 0.9 lag- 
ging power factor except for the customer, whose assumed power factor is 0.85 [4]. 
The procedures used to model the loads and source are presented with greater detail in 
Chapter IV. 
2. Field Measurements 
Harmonic levels of the system under test were measured on May 18, 1990. 
That study, hereafter referred to as the field study, measured RMS current and poten- 


tial levels of the fundamental, third, fifth, and seventh frequency components at two 


subsystem locations. Nodes two and 24 of Fig. 1 were the measurement points and V,, 
Ty 9s E5217) 19.3) 19.15, Vg, and Io4.95 were recorded for the indicated harmonics. [4] 

The voltage and current total harmonic distortion (THD) were both measured 
directly at node two; however, because of equipment limitations THD values were 
computed from the other measurements at node 24. [4] The voltage THD is defined 


to be 


(1.4) 


where V, is the RMS magnitude of the voltage fundamental and V, is the RMS magni- 


tude of the i" voltage harmonic. The current THD is defined the same way except that 
the voltage quantities are replaced by their current counterparts. [6:pp. 8] 
Measurements were recorded for six different subsystem configurations. First, 
all five subsystem capacitor banks external to the customer at nodes 24, 19, 16, 11, and 
seven of Fig. 1 were energized. Measurements were observed and recorded. Then, the 
capacitor banks were de-energized one at a time until no capacitor banks were con- 


nected. After each capacitor bank was de-energized, measurements were taken and 


recorded. The capacitor bank internal to the customer at node 25 of Fig. 1 was not de- 
energized for any of the configurations [4]. It is that set of six measurements that this 


study compares to HARMFLO simulations. 


Wi. CURRENT DISTORTION 


Power distribution systems are intended to provide a sinusoidal voltage of near 
constant magnitude for their connected loads. However, there are many types of non- 
linear loads connected to power distribution systems that draw currents which are not 
Strictly sinusoidal at the frequency of the forcing function. Such loads effectively 
inject nonsinusoidal currents into the power distribution systems to which they are 
connected. These nonsinusoidal currents are a distortion to the system current flow, 


distort the system voltage, and may be periodic or aperiodic. [1],[{3] 


A. TRANSIENTS 

Transients in the current flow of a power distribution system are aperiodic, 
discontinuous, and associated with events. These events are short in duration; the dis- 
turbances they cause in the system current flow typically decay within a few cycles of 
the system forcing function. The frequency content of transient disturbances is related 
to the natural modes of the system rather than the applied forcing function. [7] This 
Study utilizes steady state analysis and does not address transients or their effects on 


power distribution systems. 


B. HARMONICS 

Any periodic signal satisfying the Dirichlet conditions can be represented by a sum 
of orthogonal basis functions. When the basis functions are sinusoids, a Fourier Series 
results. A signal represented by its Fourier Series is a sum of sinusoids at integer mul- 
tuples of the fundamental frequency. [8] In power distribution systems, the funda- 
mental is the frequency of the applied forcing function and periodic, non-sinusoidal 


currents injected by nonlinear loads are referred to as harmonics [1]. 


1. Harmonic Sources 

There are a number of various devices that produce harmonic currents. 
Whereas all harmonic producing loads are not similar, most may be categorized as 
ferromagnetic devices, electronic power converters, or arcing devices. [3],[9:pp. 2-3] 

a. Ferromagnetic devices 

Motors and transformers are the most common power distribution system 

loads [9:pp. 2-3]. Both are essentially a coil of wire around an iron core and draw 
nonsinusoidal currents from the power distribution system. Transformers generate 
significant harmonic current levels because their cores normally saturate at the peak 
magnitude of the voltage source. Motors induce lower harmonic levels than trans- 
formers because their cores include high reluctance air gaps that do not saturate. 
Typical no-load transformer magnetizing current harmonic levels relative to the 
fundamental are shown in Table 1. [3] 


TABLE 1: NO-LOAD TRANSFORMER CURRENT HARMONICS 
Harmonic Mapnitude 


50.0% 


20.0% 





Transformer magnetizing currents are small compared to currents of other 
harmonic producing loads—typically 0.5% to 1.0% of rated load currents. Harmonic 
levels in magnetizing currents depend on voltage and loading levels. Additionally, 
currents with harmonic numbers that are multiples of three are referred to as triplen 
harmonics and cannot flow to delta connected transformers if the system to which they 


are connected 1s balanced. Therefore, triplen harmonics circulate through the windings 


of delta connected transformers in balanced three phase systems rather than returning 
through the neutral conductor. [1],[9:pp. 1, 11],[10] 
b. Electronic power converters 

Although there are many types of electronic power converters in service, 
switch mode power supplies (SMPS) and line commutated devices are being used with 
increasing regularity [3]. Personal computers, copy machines, and a host of other elec- 
tronic devices use switch mode power supplies because they offer reduced size and 
weight with increased efficiency and reliability compared to other types of power sup- 
plies. [11] 

(1) Switch Mode Power Supplies. The most common harmonic producing 


single phase load is the SMPS. Figure 2 shows a typical SMPS configuration. [11] 





Fig. 2: Atypical SMPS configuration 


The diodes conduct when the input voltage exceeds the potential across 
the capacitor C in Fig. 2. Because conduction begins near the peak of the input voltage 
waveform, the current waveform is not sinusoidal and the power factor is low. The 
harmonic content of the input ac current waveform depends on the design of the 
SMPS, the applied load, and the value of C. Of these factors, C is the most important 
because the dc ripple and harmonic levels both depend on it. [11] 

To illustrate this dependence, a typical SMPS configuration loaded 


with common electronic equipment is simulated with the SP/CE program. The entries 


listed in Table 2 depict the resulting RMS harmonic levels of the input ac current as a 


percentage of the fundamental. [11] 


TABLE 2: SMPS CURRENT HARMONICS 





(2) Line Commutated Circuits. The basic configuration of a six pulse line 
commutated circuit is shown in Fig. 3. These circuits operate as rectifiers, inverters, or 
both. Rectifiers and inverters are essentially the same, except that the polarities of E 
and v, are reversed. The current, i,, flows the same direction for both; consequently, 
the direction of average power flow is opposite for rectifiers and inverters. Steady 
state operation in the inverter mode is possible only if there is a power source on the dc 


side. [12:pp. 51-52] 
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Fig. 3: A basic commutated rectifier/inverter circuit 


The circuit of Fig. 3 is termed line commutated because there is no 
provision apart from the power distribution system to turn the thyristors off. The 
power distribution system effectively turns the thyristors off by negatively biasing 
them. The current commutates from one phase to the next when the negative bias is 
applied. [12:pp. 9-10, 587],[9:pp. 15—16],[6:pp. 20] 

Many industrial applications such as dc motor drives use six pulse 
rectifiers. Two six pulse rectifiers can be connected together to form a twelve pulse 
rectifier. [13:pp. 2-1] Twelve pulse rectifiers generate significantly smaller fifth and 
seventh harmonic components in the ac line current than six pulse rectifiers. Table 3 
lists typical harmonic levels in the ac line current as a percentage of the fundamental 


for six and 12 pulse rectifiers. [6:pp. 11-15, 22] 


TABLE 3: TYPICAL RECTIFIER CURRENT HARMONICS 
vp 







Unlike the SMPS of the previous section, the rectifier referred to here 
is a three phase circuit and does not produce triplen harmonics. Single phase rectifiers 
are used in some applications, primarily household appliances, and do produce triplen 
harmonics. Additionally, the harmonic levels depend on a host of parameters; six 
pulse rectifier parameters and their effect on harmonic current levels are discussed in 
Chapter IV. [3],[9:pp. 15-16] 

c. Arcing Devices 
Fluorescent lighting, arc welders, and arc furnaces are examples of arcing 
devices. Of these, arc furnaces cause the greatest harmonic disturbances because of 


their high power ratings. Arcing devices are essentially a voltage clamp in series with 
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a reactance connected to a single phase. Figure 4 shows a typical configuation of a 


circuit with an arcing device. [3],[9:pp. 15, 29] 


Fig. 4: A typical configuration of an arcing device circuit 


The reactance limits current flow and is normally supplied by a trans- 
former, ballast, or furnace leads. AC input current harmonic levels for typical single 
phase arcing devices are listed as a percentage of the fundamental in Table 4. [3], 


[9:pp. 15, 29] 


TABLE 4: ARCING DEVICE HARMONICS 







C. EFFECTS OF CURRENT DISTORTION 

Power distribution systems can normally absorb substantial levels of harmonic cur- 
rents without significant difficulties. Harmonic currents generally flow from their 
nonlinear load sources toward the power source. If the generation source impedance is 
small compared to the system load impedances, voltage distortion will be low and most 
harmonic currents will find a return path through the system generators. The harmonic 
currents do, however, cause increased losses in power transformers and generators. 
Additionally, localized hot spots may develop and cause insulation failure. [1],[3], 


[9:pp. 44-46] 


i 


Most power distribution systems have capacitor banks installed at strategic lo- 
cations to correct low power factors. These capacitor banks substantially change the 
nature of the system impedance. Instead of a low, purely inductive system impedance, 
the system impedance is mixed and the circuit will have one or more natural resonant 
frequencies creating low impedance resonant paths. If the resonant frequencies are 
near odd harmonic frequencies, harmonic currents may be diverted to the these low 
impedance resonant paths to neutral from their normal path through the system genera- 
tors to neutral. [1],[3],[9:pp. 19, 26-28, 43-49] 

System resonances may be series or parallel. Of the two, parallel resonances are 
more common and occur most frequently with the generation source and capacitor 
bank reactances. The combined parallel impedance at the resonant frequency is large 
and causes voltage distortion and increased harmonic current levels. Capacitor banks 
can fail from excessively high voltage and/or heat from increased current levels created 
by this phenomenon. [1],{3],[9:pp. 26-28, 43-49] 

Series resonances occur when capacitor banks are installed near feeder ends or on 
transformer secondaries. In either case, the capacitor is near the end of the line and the 
resonant impedance is the capacitor bank in series with the inductance of the line or 
transformer windings. Series resonances divert harmonic currents through the resonant 
path. As a result, loads not normally submitted to harmonics will be exposed and 
remote capacitor banks may fail. [1],[3],[9:pp. 26-28, 43-49] 

Voltage distortion is the primary problem associated with harmonics. Resistive 
loads absorb more power because of increased voltage levels. Voltage harmonics 
generate flux harmonics in motors which create losses and cause heating and 
vibrations. Additionally, circuit breakers trip unnecessarily and electronic equipment 


malfunctions. [1],[3],[9:pp. 43-49] 
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II. SOLVING THE POWER FLOW PROBLEM 


A. THE FUNDAMENTAL LOAD FLOW 

A power distribution system study based on the power flowing through the system 
is referred to as a load flow study. The active and reactive power load at each node 
and the self and mutual admittances of all the lines must be specified. These power 
and admittance inputs are used to solve for the node voltages throughout the system 
which are normally the quantities of interest. [5:pp. 193-196] 

The specified power flow to a node n is related to the voltage and current at that 


node by the relation 
P,+jQ, =V,l, (3.1) 


where j/ =V-1, P,, is the active power, Q,, is the reactive power, V,, is the voltage, and 
[,, is the current flow to node n [5:pp. 193-196]. The asterisk on / indicates complex 
conjugation. Solving for the current, 
, =a 
V 


(3.2) 


Active power is absorbed when current flows through an admittance. Apparent 
power is conserved at each node. Applying Kirchhoff's current law at each node and 
Ohm's law between each node, the voltage at any node nin an N node system is 

i 2a 
V, =| += - yy, (3.3) 
y, V, k=} 


AA 
where Y,,, the self admittance at node n, is the sum of all admittances terminating on 


node n [5:pp. 170]. The mutual admittance, Y,,, 1s the admittance between nodes n and 


k. The voltage magnitude and phase angle are normally specified at the system power 
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source called the swing bus. The swing bus is usually numbered one; therefore, V,, is 


calculated for values of n from two to N. [5:pp. 170, 193-196] 

The system of N — 1 nonlinear equations provided by equation (3.3) is of the form 
Ax = b and is frequently solved iteratively using the method of successive 
displacements, or Gauss-Seidel method. The procedure involves making an initial 
guess at the voltage of nodes two through N. These voltage estimates are inserted on 
the right hand side of the equations in (3.3) and evaluated for the voltages on the left 
hand side. The voltage calculated becomes the updated estimate. When the change 
between successive updates at all the nodes decays to a specified tolerance, the 
estimate is assumed to be the solution. [14:pp. 474 481],[5:pp. 194-196] 

If the matrix A is positive definite, convergence is guaranteed for any initial guess. 
The difference between successive updates at each node will decrease by a nearly 
constant factor, c. The number of iterations required for convergence for a tolerance € 


is at least 
m2—. (3.4) 


This number m is normally prohibitively large. Consequently, the correction to each 
update is frequently scaled by some constant called an acceleration factor. The 
calculation of an appropriate acceleration factor is difficult in the best of circumstances 
and is normally found by trial and error. For power systems, an acceleration factor of 
1.6 is frequently used [5:pp. 194-196]. [14:pp. 474-481] 

Another common way to solve the load flow equations is with the Newton- 
Raphson method of rootfinding. It is also an iterative method that requires an ininal 


estimate and is based on a Taylor series expansion of equation (3.3) in the polar form 


Ve eee) (3.5) 





P,-jQ,=>, 


k=) 
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The phase angle of Y,, is ©,, and 6, is the phase angle of the voltage at the i" node. 


[S:pp. 196-199] 
Decomposing equation (3.5) into its real and imaginary parts, the power compo- 
nents P,, and Q, are 
N 


P=) 


k=} 


V,V,Y,,|cos(O,, + 6, -5,) (3.6) 





and 


V,V,Y,,|sin(©,, {CO 5,). (3.7) 


0,=-> 


k=} 








V, 





Both P, and Q, are functions of the voltage magnitude and voltage phase angle 6.. 


Let Vv, and 6”) be the initial guess at the solution. The Taylor series expansion of 
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where the partial derivatives are obtained by differentiating equation (3.6). Evaluating 
yo 


V ,5,) for the first iteration is 


P,( 








p = Sy vVv,y,| cos(O,, +50) - 5) 4. A 
k=) 

















equation (3.6) with and 6° and subtracting it from equation (3.8), the change in 


active power for the first iteration is 


OP, 


+ Ab OF, i 
av. ' 


AP = po — pi) = eee 
1 96, |, (3.9) 


y 

















where oad ’ is the value of active power at node n that was specified when the problem 
was defined. Making the standard engineering approximation by neglecting terms of 
order two or more, the resulting system of NW — 1 equations for an arbitrary iteration i in 


matrix form is [5:pp. 196-199] 


IS 
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where the >; notation indicates >—| . 


A similar set of N — 1 equations is obtained by expanding the reactive power equa- 
tion (3.7) using the same procedure presented for the active power equation (3.6). The 


combined set of 2(N — 1) equations in matrix form is 


op : oP 
| Ap] | a5 | av] [ as® 
359 7 gv] | | (3.2 


2(N-1)x2(N-1) 
where AP", A5, and Aly are the (N-1)x1 column vectors of equation (3.10). 
The column vector AQ" is the reactive power equivalent to the active power column 
vector AP“ and the 2(N-1)x2(N—1) matrix is the jacobian, J”. The top half of 
J is the matrix of equation (3.10) and the lower half is its reactive power counterpart. 
[S:pp. 199] 

The initial guess, \v| and §), is used to evaluate equations (3.6) and (3.7) for 
P and On The differences between these calculated values and the problem 
definition values are AP and AQ. The partial derivatives of equations (3.6) and 


(3.7) are evaluated with Wee? | and §® to form J. [5:pp. 196-200],[2:pp. 12] 
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The values obtained for AP’, AQ”, and J” are used to solve equation (3.11) for 


Ad” and Aly using an appropriate numerical method. Forward backward substitu- 
tion is one common alternative to matrix inversion because it requires about one-third 


as many operations as inverting J. [14:pp. 92],[3:pp. 12] 
The values computed for Ad” and alv| are used to update the node voltages 


Naga =|V"}+alv (3.12) 





and 


58) = 6 405". (3.13) 


Updates are computed until AP and AQ decay to a specified tolerance or the maximum 
number of permissible iterations is reached. If the desired tolerance of AP and AQ is 
obtained, the final node voltages are the solution. [2:pp. 12],{5:pp. 196-200] 

The Newton-Raphson method 1s simple in form and has great flexibility. Divided 
differences are sometimes used to avoid computing partial derivatives. The number of 
iterations for convergence depends on the initial guess. If the initial guess is far from 
the solution, the iterate may diverge. Conversely, when the iterate is sufficiently close 
to the actual solution, convergence is rapid. Each iteration requires at least N?+N 
function evaluations for a system of N nonlinear equations making it a very 
computationally expensive method. For that reason, other methods such as Gauss- 
Seidel may be used to refine the initial guess before applying the Newton-Raphson 


method. [14:pp. 54-56, 92-95] 


B. THE HARMFLO SOLUTION 
1. Solution Method 
The HARMFLO problem formulation differs from the fundamental load flow 


problem of the previous section because it includes harmonic generating loads. 
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Detailed load models for line commutated circuits and nonlinear resistors are included. 
The terminal voltages of the defined nonlinear loads are represented by their Fourier 
Series and are related to their input currents by these models. [2:pp. 12-34] 

The nonlinear loads are treated like current sources connected to the power 
distribution system. The initial estimate of currents injected and terminal voltages for 
nonlinear loads are obtained from model data based on the specified loads. Because 
reactances are frequency dependent, the power distribution system has self and mutual 
admittance matrices for the fundamental and also for each harmonic frequency of the 
injected current. Although these admittance matrices only relate the currents at a given 
frequency to the voltage at that same frequency, the system voltage response to one 
harmonic frequency may be related to the current injected at a different frequency. 
Additionally, the voltage distortion caused by injected harmonic currents will affect the 
injected current making an iterative solution necessary. [2:pp. 12—34],[17:pp. 270-283] 

The iterative Newton-Raphson method is reformulated to include harmonic 
frequencies. Additional equations are required to solve for the harmonic node 
voltages. The equations are based on Kirchhoff's current law and the conservation of 
apparent power where appropniate. Xia and Heydt develop the additional equations 
required to accommodate line commutated circuits in Refs. [15] and [16]. The 
HARMFLO formulation extends these equations to include nonlinear resistors and a 
variety of parameter combinations for line commutated circuits not included in the 
original development. [2:pp. 12-34] 

Whereas detailed formulation of the equations is beyond the scope of this 
study, the general results are stated here for completeness. For an n node system with 


m nonlinear busses, the linear busses are numbered one through m-— 1. The swing bus 
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is the linear bus numbered one. The nonlinear busses are numbered m through n. 
Non-triplen odd harmonics from five to L are considered. [2:pp. 11-16] 


The active and reactive power balance is 


[aw]=[J5 lee J }fav® AVS)... avy (3.14) 


ere 
AQ 


order. The fundamental jacobian, J“, is the same as the square matrix in equation 


Ad 
| [av | i and the superscript k indicates harmonic 


where [AW]= 
(3.11). Partial derivatives of equations (3.6) and (3.7) evaluated with ” harmonic 
frequency component values are used to build the harmonic jacobian J™. [2:pp. 11-16] 

The nonlinear device models use two state variables, @ and B. The change in 


state variables for an iteration is defined to be 


[A®] =[ Aa, AO ns a Aa, AS, AB 1 oie As, ]'.(3.15) 


(&) 


[7 


The k” harmonic current injected at node f has real and imaginary parts 2’ and ae 
respectively, where St<n. The partial derivatives of nonlinear device currents with 


respect to nonlinear device state variables at the k” harmonic are [2:pp. 11-16] 


ag") ag) 
. | oa, of, 
da, op, 


The harmonic jacobian that relates the K” and j” harmonics is indicated by 


YG") and is defined as 


(4,4) (kk) pj 
yGrrasJ¥  t@ kay (3.17) 
Gi) ke j 
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where Y“*) is a matrix containing partial derivatives of the kK” harmonic injection 
currents with respect to the k” harmonic bus voltages derived from the system 
admittance matrix. The partial derivatives of the k* harmonic device currents with 
respect to the /“ harmonic applied voltages are derived from the nonlinear device 


models and form the matrix 


O 46n-1),2(m-1) 2(m-1),2A 
need : Ogry ag’) (3.18) 
O2,.2(m-1) - VA GQ" ) ov) 
(On (z) (k) 
og, f og, t 


Vg) yA 


where V,“) and © are the K” harmonic voltage magnitude and phase angle at the f” 
bus, O; , is an ix j matrix of zeros, and mSt <n. [2:pp. 11-16] 
If h harmonics in addition to the fundamental are considered, the set of 2n(1 + 


h) + 3m nonlinear equations in matrix form is 


AW Jy” J nee J (L) 0 AV”) 

AY YG YE (1,5) — YG) H” AV’) 

AY) |=| YES" yGS) =. v¥GCAE Ae on ee 
ee : AV“) 

| AI” | | yGe» yqe» = YG) yo | A® 


These equations are solved in much the same way as the fundamental load flow equa- 
tions of the previous section are solved. An initial guess is made at the harmonic 
voltages and nonlinear device state variables. The power and currents are evaluated 


with the estimated values. Changes in voltages and state variables are computed from 
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equation (3.19). With each iteration, the voltage and state variable estimates are up- 
dated with the computed values until the change in power and currents decay to the 
specified tolerance. [2:pp. 11—16] 

The solution speed and memory requirements depend heavily on the system 
because of the large number of sparsely populated matrices. Typical solution times for 


two system sizes are shown in Table 5 for three common machines. 


TABLE 5: TYPICAL SOLUTION SPEED FOR THREE COMMON MACHINES 


aes speed 
CDC 6600 | Cyber205 | VAX 11/780 
69 bus, 7 harmonics 800 seconds 3200 seconds 
21 bus, 7 harmonics 350 seconds 1400 seconds 






A HARMFLO version is also available for the IBM PC. The PC version is pre-compiled, 
dimensioned for a 50 bus system, and requires 512 kilobytes of main memory and a 
math coprocessor. A 100 bus version can be compiled in about two hours and requires 
an estimated one megabytes of expanded memory. [3:pp. 3.1—3.3, 3.36—3.38] 

2. Assumptions and Limitations 

In addition to the 50 bus limit for the PC version, the HARMFLO problem for- 
mulation assumes a balanced bilateral three phase system and delta connected nonlin- 
ear devices. Consequently, single phase loads are not permitted and all even and 
triplen harmonics are excluded. [2:pp. 12-13] 

Nonlinear loads may only be one of the defined types. A static transformer tap 
ratio or phase shift may be defined; however, tap or phase changing under load is not 
allowed. The frequency range is limited to 1860 Hz and the fundamental can only be 
60 Hz. Additionally, if the apparent power is specified at a nonlinear bus instead of 
the reactive power, the reactive power is computed without any range limits. [13:pp. 


3.1, A.1-A.7, 1.6] 
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IV. MODELING THE SYSTEM UNDER TEST 


A. THE POWER SOURCE 

A 69 kV substation bus draws power from a station bus. The substation bus 
potential is stepped down to 12.5 kV and distributes power to Audrain County. The 
Thevenin equivalent impedance at the Audrain 12.5 kV feeder is 0.947 per unit. This 
figure includes the short circuit impedance of the transmission system and the sub- 
transient synchronous reactance of the system generators. A one line drawing of the 


system power source is shown in Fig. 5. [4] 


0.94} p.u. Cs 






Audrain County Feeder 






Station Bus Substation Bus 


12.5 KVAC 
@) 


Fig. 5: A one line drawing of the system source 


With HARMFLO, a system is driven by a source called a fundamental swing bus. 
The bus voltage is controlled to a percentage of the system voltage base. Only funda- 
mental frequency currents are allowed on the swing bus. Once the source is specified, 
it may be connected to the rest of the system. [13:pp. 3.1—3.12] 

Connections are specified by which nodes they connect and the equivalent 
impedance of the line expressed as a percentage of the system impedance base. This 
study does not address the skin or long line effects; however, HARMFLO does have 
hyperbolic long line and resistive skin effect models. These options are invoked by 


entering the base impedance and line length parameters. [13:pp. 3.14—3.18] 
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In addition to the normal paths, an option is provided for a harmonics only path. 
This option is necessary because harmonic currents will not have a path to ground 
through the swing bus. A path shunt to the swing bus with impedance equal to the 
short circuit impedance plus the subtransient synchronous reactance is required. The 
HARMFLO equivalent circuit model for the system source is shown in Fig. 6. [13:pp. 


3.5-3.7] 







j94% 


100 / 0° % 


60 Hz Path 2 to 0 for 


onics Only 


= © 


Fig. 6: The HARMFLO model of the system source 


Note that the harmonic path to ground is connected to node zero. The connection be- 


tween node zero and ground is implicit. 


B. THE SOURCE OF HARMONICS 

The dominant sources of harmonics in this system are the seven 100 HP adjustable 
speed dc drives shown lumped together at node 25 in Fig. 1. DC drives require ac to 
dc converters and rectifiers are the principle elements of such converters. Before simu- 
lating the entire Audrain subsystem, a circuit model of the dc drive and system source 
is simulated to determine the dc drive parameters. Parameters are adjusted until the 


desired response is obtained. 
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1. Modeling the Rectifier with HARMFLO 
Six pulse rectifiers are used for this study because they are common to indus- 


trial applications such as dc motor drives. The HARMFLO circuit model of a six pulse 


rectifier is shown in Fig. 7. [13:pp. 2.1—-2.4] 


PTET MPSS ed eeseceens Soe ee ee SOeBeese Bee saree ony 





Fig. 7: The HARMFLO circuit model of a 6 pulse rectifier 


Rectifiers may or may not have a dc source on the dc side of the circuit in Fig. 
7. For dc drives, this source is equivalent to the back electromotive force voltage of 
the motor. If the dc source is chosen, HARMFLO will calculate the dc speed voltage E 
as one of the device state variables. If a rectifier without a dc source is chosen, E = 0. 
Frequently, a dc filter capacitor is used in shunt with the dc load. However, the effect 
on the harmonic content of the ac line current is small and the HARMFLO model does 
not include it. The parameters X, R, F, P, and Q identified in Fig. 7 must be specified 
as input parameters to HARMFLO. [13:pp. 2.1—-2.4, B1-B10] 

Rectifier transformers are modeled with HARMFLO by the equivalent reactance 
of the windings, X, expressed as a percentage of the system power base. The connec- 
tion type such as Y~Y or Y—A must also be specified. Although the transformers are in 


three phase circuits and are not allowed to pass triplen harmonics with either connec- 
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tion type, the voltage phase angles differ by the standard 30° between Y and A connec- 
tions. Additionally, if the resistance of the windings is significant, it may be lumped 
with the other components that make up the equivalent resistance, R. [13:pp. B1-B10] 

For dc drives, R is comprised of the rectifier circuit and armature resistances. 
The reactance F represents the equivalent circuit reactance including filters and the dc 
machine reactance. Both are expressed as a percentage of the dc impedance base. The 
fundamental active and reactive power load of the drive are represented by P and Q, 
respectively. A permissible alternative is to specify the apparent and active power. In 
either case, the power quantities are expressed as a percentage of the system power 
base. [13:pp. 3.1-3.33] 

2. DC Drive Parameters 

The field study indicates that the customer operates the seven 100 HP 
adjustable speed dc drives on a one MVA base with a load of 500 kW and 309 kVAR 
[4]. Initially, one 700 HP adjustable speed dc drive is considered. One drive requires 
one six pulse rectifier for a dc voltage source. The rectifier parameters R, F, and X are 
unknown. Several assumptions are made to estimate these parameters for input to 
HARMFLO. 

Assuming the rectifier transformer steps the 12.5 kV supplied by the utility 
down to 440 V and a commutation angle a = 0, the average dc voltage is 


V, a7 === (440 V) = 594.21 V. (4.1) 


342 3/2 
=——-V, cos 
T 
For a 700 HP load, the dc current is [12:pp. 46] 


|, =—=———-_ = 841.46A. a 
“Vv 594.21V Gua 


25 


Assuming the converter has a dc filter choke with 5% losses, the loss due to the 


choke is 


I?R=0.05x 500 kW =25 kW. 
Therefore, the choke has resistance 


Poe 25KW 


= mQ. 
I? (841.46 A) 


(4.3) 


(4.4) 


Since the type of wire in the choke is unknown, number 4/0 braided copper 


wire is assumed because of its current rating. Number 4/0 copper wire braided 19 by 


0.106 inches has resistivity 0.049 ohms per thousand feet, a 0.528 inch copper diame- 


ter, and a 0.1662 square inch copper cross section. Polyurethane or polyvinylchlonde 


insulation adds 550 pounds to every 1000 feet of wire for every cross sectional square 


inch of conductor. [18:pp. 277] 


The density and cross sectional area of the insulation are 





(so) 
100 ft 
Pins =X 0.1662 in’ = 0.09141 2s 
in ft 
and 
Ps 0.091428 
3 =) [Ay = | eee: In =U. in. 
A, =| FIA = 0.1662 in? = 0.0237 in? 
Pre 0.6405—— 


Therefore, a number 4/0 insulated wire has a diameter of about 


(4.5) 


(4.6) 


a) : 2 
dyay = 2 Ain a 2 =) [Soest (ee) =0.556in. (4.7) 
y It 
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The choke has a total wire length of 


WoL = us ae = 720.61 ft = 219.63 meters. 


p ee . (4.8) 





Assuming a choke diameter of one meter, the circumference is 3.14 meters. With the 


assumed diameter, the choke has 





W.L. 219.63 
= ——___——_-= = 70 turns (4.9) 
circumference N 
and a corresponding physical height of 
h =(70)(0.556 in) = 0,99 meters. (4.10) 
in 


The inductance of an air core solenoid choke is 


_ HyN*A 


ie 4.11 
h (4.11) 
This expression can be manipulated to be 
W. L;. 2 
mad) Hol Se) _ pat tr (4.12) 
Nd d Id 


where r is the radius of the solenoid and d is the diameter of the wire.[19:pp. 268] 


Numerically, the choke has inductance 


L=—e = = 4.9 maz. (4.13) 


ay 


The dimensions and inductance are reasonable for an air core solenoid choke used to 
filter the converter of a 700 HP dc drive. Figure 8 illustrates the assumed dimensions 


of the choke and the wire used to construct it. 


#4/6 Cepper Wire 4.9 mH Air Cere Seleneid Cheke 
Braided 19 X 6.166 inches N = 70 Turns 
with PE er PVC insulation < Cc \ 





Cepper Area 0.1662 in : 
Dimensions in Inches Dimensiens in Meters 


Fig. 8: Wire and choke dimensions 


The motor base impedance is 


4 2 
gt = Oo ee (4.14) 
S S8&8kVA 


Assuming the armature resistance of the dc drive is 5% of the motor base impedance, 


the armature resistance is 


R , = 0.05 x329.3 mQ = 16.5 mQ. (4.15) 


The forward drop resistance of the rectifier SCR thyristors is neglected; therefore, the 


total resistance on the dc side of the circuit in Fig. 7 is 


Ry, 


0 


wat = Ringe +R, =35.3mQ + 16.5 mQ = 52 mk, (4.16) 
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To convert the rectifier parameters for input to HARMFLO, the rectifier trans- 
former secondary is assumed to be 440 V. The base dc voltage is 


Vane = Veacu-n = a = 254.03 V. (4.17) 


For consistency, the dc impedance base is defined to be the same as the ac impedance 


base and is given by 


eh) Se = 193 6 nO. (4.18) 


Note that the appropriate power base is the one MVA power base which the customer 
operates on and not the 100 MVA system power base. [13:pp. B6—-B9} 
The equivalent resistance, R, is 


R 52 m2 
Pepe eta) = 2019686 pu. 4. 
Eva 7 1OsGunO 3 ow) 


This per unit value must be further converted to the percentage 
R, = 100 R,, = 26.86%. (4.20) 


The motor reactance is small relative to the choke and is neglected. Consequently, the 


total reactance in the dc circuit is [13:pp. B6-B9}] 





X, =120aL =1207(4.9 mH) = 1.84 Q, (4.21) 
and 
pate oA 109 = 951.62 %. (4.22) 
Zac 0.19362 


The rectifier transformer reactance is assumed to be 10% of the customer's 
impedance base. Several combinations of Y and A connections are evaluated for 
model performance. The transformer reactance expressed as a percentage of the 


system base is 
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6 
oy X10 07 ON (4.23) 
100 x10 
The given rectifier load must be converted to a percentage of the system power 
base. The active power load is 


_ 500 KW 


A = OMe 50 %. (4.24) 


The reactive power load obtained using the power triangle is 


[o2_ p2 588,000? — 500, 0002 
VST =P 100 = ¥288:000" = 500, 000" 100 =0.31%, (4.25) 


Os =e 100 x10° 


With these source and rectifier parameters, the minimal circuit model of Fig. 9 
is input to HARMFLO and simulated. Model parameters are then varied and model 


performance assessed. 
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Fig. 9: Minimal single phase circuit model 


3. Rectifier Model Performance 
a. Theoretical Behavior 
The rectifier will generate harmonics in the ac line current depending on the 
delay angle, @, and the commutating reactance, X [6:pp. 22]. The line current drawn 


by a six pulse rectifier is typically a quasi square wave. Figure 10 shows the theoreti- 
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cal line current for a six pulse rectifier with a Y-Y transformer connection. [12:pp. 57— 


59] 


oe | ae Line Voltage (V, ,) L 
é Galesie os > SS ° | 


é ad: / 
f Lewes Z 
é Q) : Line Current (I,) f 
7 * 3 
*, Fs 





Y - Y Transformer Connection  *., 


0 0.01 0.02 


Time (seconds) 
Fig. 10: Theoretical six pulse rectifier ac line current 


The line and phase voltages of a Y-Y connected load are in step. Thirty 
degrees separates the negative to positive zero crossing of a sinusoidal phase voltage 
from the crossover point with an adjacent phase voltage. The angle at which the start 
of the quasi square wave exceeds this 30° is defined as the delay angle, a. [12:pp. 55— 
60] Whereas the commutating reactance, X, must be specified by the user, @ is nor- 
mally computed by HARMFLO as one of the nonlinear device state variables. However, 
the delay angle may be specified as an input if it is known. [13:pp. 3.28-3.33] 

The rising and falling edges of the quasi square wave indicated by bullet 
one in Fig. 10 are a function of the rectifier transformer reactance, X. If X = 0, the 


slope of the rising and falling edges is infinite and the commutation interval is zero. If 
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X is non-zero, the commutation interval, u, 1s non-zero. The commutation interval is 
also referred to as the commutating angle, U. A non-zero commutating angle produces 
an observable notch in the line voltage because two adjacent phase voltages are shorted 
together through an SCR and the commutating reactance of the rectifier transformer 
during each time interval u. For this reason, the commutation interval, u, is also 
known as the notch width for the line voltage. [12:pp. 55—60],[6:pp. 11-15] 

The flat surfaces of the quasi square wave indicated by bullet two in Fig. 10 
are a function of the dc reactance, F. If F is infinite, the flat surfaces have no ripple. 
For finite values of F, the nipple depends on the value of F. The value of R determines 
which values of F satisfy the convergence criterion. [12:pp. 55—60],[6:pp. 11-15] 

b. Actual Performance 
Rectifier line current components are obtained from HARMFLO output for 
the circuit of Fig. 9. A time domain signal of the form 
i(t) = > M, sin(120ant+ ¢,) (4.26) 
n=] 
is constructed where M, is the magnitude and 9, is the phase of the n” harmonic 
current. The harmonic number nv is considered for the values n = 1, 5, 7, 11, and 13. 
Simulated rectifier line currents for the circuit of Fig. 9 are shown in Fig. 11. 

Model performance is identical for Y-Y and A—A connected rectifier trans- 
formers. It is also the same for Y-A and A-Y connections. Additionally, the mixed 
transformer connection models all produce identical magnitudes and phase angles to 
the like connection models with two exceptions. The phase angles of the 5" and 7% 
harmonics are shifted by 180°. This difference is the result of the 30° phase shift be- 
tween phase voltages of Y and A connected transformers. Magnitudes of the current 


harmonics relative to the fundamental are listed in Table 6. 
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TABLE 6: SIMULATED HARMONIC LINE CURRENT MAGNITUDES FOR THE 
CIRCUIT OF FIG. 9 (RELATIVE TO THE FUNDAMENTAL) 


Tiger | sitet | see 
Harmonic # with a de source without a dc source from Table 3 






The simulated values of Table 6 are considerably different than the typical 
values of Table 3. To isolate the cause, each parameter is varied independently. The 
rectifier transformer reactance X is varied over a range from 0.01% to 1000%. The 
simulated results for the model with a dc source more closely match the typical values 
as X is increased. With X = 100%, the fifth harmonic is only 4.4% smaller relative to 
the fundamental, whereas at 1000%, it is 23.2% less. The model without a dc source 
shows little variation as the transformer reactance 1s changed. The fifth harmonic cur- 
rent magnitude changes only 1.4% relative to the fundamental over the full range of X. 


Figure 12 shows the extreme case of X = 1000% with Y-Y connected transformers. 
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2 
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Fig. 12: Simulated rectifier line current with extreme transformer reactance 
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The dc circuit reactance, F, is varied over a wide range. The rectifier line 
current for the model without a dc source does not change; however, the line current 
for the model with a dc source has less mpple when F is increased. The HARMFLO 
solution converges for values of F between 100% and 1350% with R = 26.86%. Fig- 


ure 13 illustrates how the line current changes for the model with a dc source. 


1.5 


Current (per unit) 





F = 100% (0.5 mH) ==> { 


0 0.01 0.02 
Time (seconds) 


Fig. 13: Simulated rectifier line current as a function of de circuit reactance 


The equivalent dc circuit resistance, R, 1s varied over a wide range and 
makes no significant difference in the line current for the rectifier model with a dc 
source. For the model without a dc source, R is calculated by HARMFLO from the 
voltage and active power at the rectifier bus. The active power load level affects the 
line current as illustrated in Fig. 14. 

The power factor also affects the rectifier line current. Figure 15 shows 


how the line current changes in response to a varying power factor. 
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(b) Rectifier without a de source 


Rectifier line current as a function of active power load 


Fig. 14 
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(a) Rectifier with a de source 
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(b) Rectifier without a de source 


Fig. 15: Rectifier line current as a function of power factor 
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The model that best approximates the quasi square wave of Fig. 10 is the 
rectifier with a Y-Y or A—A transformer connection without a dc source. Simulation 
results indicate that the estimated parameter values are reasonable. However, the fifth 
harmonic of the line current is significantly larger than typical values. Furthermore, 
the customer operates seven drives. It is likely that the drives are operated indepen- 
dently. Therefore, modeling the drives by lumping them all together is an inadequate 
approximation. 

4. A Revised Rectifier Model 

For a system with multiple nonlinear loads, harmonic levels may be reduced by 
cancellation. This cancellation is achieved by alternating between Y—Y and Y-—A rec- 
tifier transformer connections at adjacent loads. If the loads are split evenly between 
Y-Y and Y—A connected converters, a 12 pulse rectifier results and the fifth and sev- 
enth harmonic line currents are significantly reduced. [13:pp. 2.8] 

The seven 100 HP adjustable speed dc drives are split between Y-Y and Y—-A 
connected rectifiers. However, there are an odd number of drives making it impossible 
to split the load evenly. Assuming the customer powers four of the drives with Y—A 
connected rectifiers and three drives with Y-Y connected rectifiers and that the load is 
evenly distributed among all seven drives, a 57% Y—A to 43% Y-Y load split results. 
The parameters X, R, F, P, and Q are all scaled accordingly. The revised minimal 
equivalent single phase circuit model is shown in Fig. 16. 

The circuit of Fig. 16 is simulated for rectifiers with and without a de source. 
AC line currents are generated with components obtained from HARMFLO output and 
equation (4.26). The resulting waveforms are shown in Fig. 17. Harmonic levels ex- 


pressed as a percentage of the fundamental are listed in Table 7. 
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Fig. 16: Revised minimal single phase circuit model 
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Fig. 17: Simulated rectifier line current for the circuit of Fig. 16 
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TABLE 7: SIMULATED HARMONIC LINE CURRENT MAGNITUDES FOR THE 
CIRCUIT OF FIG. 16(RELATIVE TO THE FUNDAMENTAL) 


Rectifier Rectifier Typical 6 Typical 12 
Harmonic # with a de without adc | pulse values | pulse values 
source source from Table 3 | from Table 3 


| Sth 9.2% 46% | 












The simulated current waveforms shown in Fig. 17 are similar to idealized 12 
pulse rectifier line currents [12:pp. 387-389]. The simulated harmonic levels listed in 
Table 7 fall between the typical values for six and 12 pulse rectifiers from Table 3. 
The model is significantly less accurate at frequencies above the seventh harmonic. 
This inaccuracy is of little consequence for this study because field measurements were 
not recorded for harmonic orders larger than seven. 

The rectifier with a dc source is selected because it more closely resembles the 
physical circuit being modeled and also because of its fifth harmonic line current 
magnitude. At 9.2%, it is roughly half of the typical six pulse value listed in Table 3. 
A realistic goal when alternating rectifier transformer connection types might be to re- 


duce the largest single harmonic current magnitude by a factor of two. 


C. LOAD MODELING 

The computer analysis uses the estimated load values from the field study. The 
general procedure is stated here for completeness. The customer's nonlinear load is 
based on a 700 HP load with a 0.85 power factor. The difference between the total 
measured voltampere load at node 25 and the customer's estimated nonlinear load is 
the linear load estimate at node 25. [4] 

To estimate the other loads, each feeder is considered separately. The ratio of the 


total measured apparent power to the sum of all the connected transformer ratings is 
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computed for each feeder. This ratio scales the transformer ratings and produces an 
estimate of the loads connected to that feeder. For feeder 53, the scale factor is com- 
puted without the customer's load because the load at node 25 is already estimated. [4] 
Several nodes in Fig. 1 represent more than one transformer. Although there are 
some residential loads on feeder 51, feeder 52 is predominantly residential loads. A 
number of the nodes on feeders 51 and 52 represent small residential areas with more 
than one house connected. Whereas each house normally has its own transformer, the 


system grows arbitrarily large if the nodes are not grouped in this manner. [4] 


4} 


V. SIMULATION RESULTS 


The system of Fig. 1 is simulated using the revised rectifier model to approximate 
the customer's seven 100 HP adjustable speed dc drives. Six system impedance con- 
figurations are simulated. These configurations facilitate evaluation of nonlinear load 
modeling accuracy and prediction accuracy of system-wide harmonic propagation. 


Table 8 lists the six system configurations. 


TABLE 8: SIMULATION IMPEDANCE CONFIGURATIONS 


7 


Capacitors 
Connected 
at Nodes 





The accuracy of the simulated current at node 24 is a measure of the nonlinear load 


model effectiveness. Table 9 lists the simulation results and field measurements of /,,. 


TABLE 9: CURRENT AT NODE 24—FIELD DATA AND SIMULATION RESULTS 
(Amperes) (% of Fund.) (% of Fund.) (% of Fund.) 
Data lation Data lation Data lation Data lation 
| 2 | 6892 | 55.43 | 623 | 593) | 34) 66 | ese een 
| 3 | 67.39 | 55.89 | 6.510) 5.59" (1302s eter cree 
| 4 | 6892 | 55.89 | 637 | 5.45 | 114 | 265 | 649 | 6.04 | 


| 5 | 6143 | 56.35 | 8.05 | 5.33 | 256 | 264 | 845 | 5.96 | 
|_6 | 67.39 | 5681 | 6.10 | 5.18 | 2339) 250 | eee eae 







The simulated current fundamentals are consistently smaller than the measured values. 


However, the simulated harmonic current to fundamental current ratios are similar to 


42 


the measured ratios. Figure 18 illustrates how the simulated and measured values 


change in response to a varying system impedance. 





(a) Field measurements 


ese 230 9 
Go esp VY WO FSF tA DW ~J CO 1 





(b) Simulated values 


Fig 18: Node 24 current harmonics and THD relative to the fundamental 
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The simulated and measured voltage quantities at node 24 are similar. Table 10 


summarizes the simulated and measured values. 


TABLE 10: VOLTAGE AT NODE 24—FIELD DATA AND SIMULATION RESULTS 


(kV) (% of Fund.) (% of Fund.) (% of Fund.) 
Data lation Data lation Data lation Data lation 
| 1 {| 1240 | 1257 | 096 | 160 | 060 | 034 | 113 | 1.64 _ 
| 2 | 1230 | 1245 | 034 | 097 | 0908 | 058) 01975| ee 
| 4 |{ 1230 | 1229 | 031 | 9059 | 0.75 | (oe) esi) oes 
| 6 | 1230 | 1218 | 033 | 048 | 027 | 032 | 043 | 0.58 _ 







Unlike the current fundamentals at node 24, the simulated voltage fundamentals are 
not significantly different than the measured values. Figure 19 illustrates the differ- 


ences in simulated and measured voltage THD. 


Field Data 
Simulation Results 





Run 


Fig. 19: Simulated and measured voltage THD at node 24 
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Comparisons of the simulated feeder currents to field measurements indicate how 
closely the system model approximates the actual system and how well HARMFLO pre- 
dicts propagation of the injected harmonic currents. Tables 11, 12, and 13 summarize 


the simulated and measured current values for feeders 51 through 53. 


TABLE 11: FEEDER 51 CURRENT, ,,—FIELD DATA AND SIMULATION RESULTS 
(Amperes) (% of Fund.) (% of Fund.) (% of Fund.) 
Data lation Data lation Data lation Data lation 
| 1 | 23.74 | 2494 | 136 | 330 | 117 | 1.18 | 186 | 3.51 
ee e554 e4e4 145 | 97 | 217 | 200 | 254 | 281 
eo | 268772494 | 1401.35) 245 | 252 | 282 | 284 | 
m4 | 2637 | 25.40 | 158 | 115 | 193 { 197 | 250 | 230 | 


| 5 | 29.79 | 27.71 
6] 31.13 | 27.71 | 019 | 0.26 | 030 | 0.18 | 120 | 0.32 












Fundamental Fifth Seventh 
(Amperes) (% of Fund.) (% of Fund.) (% of Fund.) 


Data lation Data lation Data lation Data lation 










Data lation Data lation Data lation Data lation 
1_ | 134.97 | 109.93 | 5.20 | 7.64 | 147 | 0.70 | 640 | 7.66 | 


| 2 | 143.74 | 108.08 | 3.34 | 465 | 1.29 | 1.21 | 450 | 481 | 
| 3 | 146.15 | 113.63 | 3383.08 (060 | 144 | 430 | 3.36 _ 
, | 1.04 


(Amperes) (% of Fund.) (% of Fund.) (% of Fund.) 
|. 







79.86 | 2.50 | 1.04 
_ 5 | 151.73 | 119.17 | 283 | 247 | 1.00 | 1.15 | 380 | 272 | 
_ 6 | 151.73 | 120.09 | 259 | 241 | 097 | 116 | 350 | 268 | 
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Figure 20 shows how the simulated and measured THD values compare for the 


three feeders. 
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to the fundamental for feeders 51 through 53 


Current THD relative 


Fig. 20 
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Referring to Tables 11-13, the simulated feeder 53 current fundamental is signifi- 
cantly smaller than the field measurements for all system configurations except run 
number four. The simulated feeder 51 and 52 current fundamentals are close to the 
measured values. Additionally, the simulated and measured harmonic currents and 
current THD are within a few percent in most cases and vary the most on feeder 52. 

Perhaps the best overall indicator of the system modeling effectiveness is the com- 
parison of measured and simulated voltage at the Audrain County 12.5 kV feeder trunk 
and the current supplied to the substation. Tables 14 and 15 summarize the voltage 


and current quantities at node two. 


TABLE 14: VOLTAGE AT NODE 2—FIELD DATA AND SIMULATION RESULTS 
el (kV) _ (Zo of Fund. ) (Zo of Fund. ) (Zo of Fund. ) 
Data lation Data lation Data lation Data lation 
| 1 | 1257 | 1260 | 050 [| 113 [ 040 | 029 | 084 | 1.16 | 
[2 | 12.47 [| 1253 | 032 | o68 | 060 | 049 | 0.81 | 084 
me 24 ie4emi 023 | 0.47 | 062 | 062 | 092 | 078 | 
| 4 | 1243 | 1242 | 021 | 040 | 059 | 050 | 0.78 | 0.64 | 


5 | 1245 | 12.38 | 020 | 037 | 040 | 039 | 062 | 0.54 
6 | 12.44 | 1231 | 017 | 030 | 0.18 | 020 | 056 | 0.36 


(Am i eres) (Zo of Fund. ) (Zo of Fund, ) (Zo of Fund. ) 
| 


















Data lation Data —- Data lation Data lation 


1__| 175.75 | 161.20 
| 2 | 180.56 | 157.05 [4.26 | 3.10 | 2.21 | 4.60 | 4.80 | 







Tae as | | 

3 [1es.55 | 15843 | 27a | 209 [280 {275 [440 {3.99 
251} 246 _ 
228 


ae 191.73 
[5 | 196.54 | 164.44 | P2171 [094 [65 {320 [2.73 


| 6 | 210.90 | 174.60 | 160 | 165 | 0.07 | 0.79 | 230 | 183 | 


Figure 21 compares the simulated and measured current and voltage THD at node two, 


the Audrain County 12.5 kV feeder. 
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(b) Voltage THD at node 2 
lated and measured current and voltage THD at node 2 


imu 
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Fig. 21 


VI. MILITARY STANDARDS 


The interface standard for shipboard systems, MIL-STD—1399 (NAVY) SECTION 
300A of 13 October 1987, entitled Electrical Power, Alternating Current (Metric) pro- 
vides guidance about power quality. Because of the assumptions and limitations of the 
HARMFLO program, comparisons between simulation results and the guidance provided 
by this standard to 400 Hz systems, current and voltage imbalances, and transient phe- 
nomenon are inappropriate. 

Specific guidance that does apply to this study includes 5% maximum voltage 
THD, 3% maximum single voltage harmonic, and 5% maximum voltage deviation 
factor [20:pp. 8]. The voltage deviation factor is defined to be the maximum differ- 
ence between the voltage waveform and a pure 60 Hz sine wave with amplitude equal 
to the maximum waveform value, relative to that maximum value [20:pp. 5]. 

Equipment must not induce single harmonic line currents greater than 3% of the 
full load rated fundamental for harmonic order between two and 32. If the power 
source characteristics are not specifically known, pulsed loads may not exceed 70 
kVA. If the source characteristics are known, pulsed loads up to 25% of the generator 
rating may be used depending on the power factor. [20:pp. 15, 35] 

Shipboard 60 Hz power distribution systems are designed to operate with power 
factors between 0.80 and 0.95. Although line-to-neutral filter capacitors are allowed, 
the preferred capacitor bank installation is line-to-line. In any case, installed capacitor 
banks may not exceed 0.1 pF for 60 Hz systems. Additionally, equipment rated at 5 
kVA or more must use a three phase 440 V input and source impedances are assumed 


to be 5 22 at 100-200 kHz unless specified otherwise. [20:pp. 14] 
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VII. CONCLUSIONS 


Relatively little is known about the parameters used to describe the system of Fig. 
1. However, most of the simulation results are remarkably close to the field measure- 
ments. The major disparities are the differences between simulated and measured cur- 
rent fundamental drawn by the customer and the current harmonics on feeder 52. 

There are at least three possible reasons why the rectifier model for the seven 100 
HP adjustable speed dc drives performs adequately with the test case and produces 
such a small current fundamental when connected to the entire system. A more de- 
tailed system model may be required to obtain the proper interaction between compo- 
nents. Another possibility is that the estimates of the customer's loads are inaccurate. 
Finally, the rectifier model internal to HARMFLO may not perform properly. Of the 
three, 1t is most likely that the system model requires more detail. 

The disparity in current harmonics on feeder 52 is probably due to the fact that 
many household appliances produce harmonic currents. The system model does not 
include any nonlinear loads on feeder 52. Additionally, grouping residential loads to 
limit the system model size may be at the expense of accuracy. 

The THD plots for the six system configurations show that there are system reso- 
nances and that HARMFLO simulations do predict system resonances with reasonable 
accuracy. The trends exhibited by the simulated and measured values are similar. 
Long line modeling and inclusion of some nonlinear loads on feeder 52 would likely 
improve results in this area. 

From the military perspective, HARMFLO has apparent applications where large 


nonlinear loads are planned for shipboard installation. The design process for large 
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rectified loads would be enhanced by HARMFLO simulations, particularly since the pa- 
rameters for such systems are more precisely known than the system of this study. 
Reference [20] should be updated to give more specific guidance about large rectified 
loads. 

The Department of Defense (DoD) has chosen ADA as the preferred computer 
language for programs used by the DoD. Potential military users of HARMFLO should 
be aware that it is written in FORTRAN. Currently, there are no programs written in 


ADA with capabilities similar to HARMFLO. 
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0001 


APPENDIX 
INPUT DATA SET FOR THE SYSTEM UNDER TEST 


AUDRAIN COUNTY MISSOURI. FEEDERS 51 - 53. 


0002 


SIMULATION OF AUDRAIN COUNTY DISTRIBUTION SYSTEM FOR 


0002 


COMPARISON TO FIELD DATA. 


0002 


ALL CAPACITORS CONNECTED 


2003 
1 swing busO1 


2 audrain 12kv03 


3 feeder 52-103 
4 feeder 52-203 
5 feeder 52-303 
6 feeder 52-403 
7 feeder 52-503 
8 feeder 52-603 
9 feeder 52-703 
10 feeder 52-803 
11 feeder 51-103 
12 feeder 51-203 
13 feeder 51-303 
14 feeder 51-403 
15 feeder 53-103 
16 feeder 53-203 
17 feeder 53-303 
18 feeder 53-403 
19 feeder 53-503 
20 feeder 53-603 
21 feeder 53-703 
22 feeder 53-803 
23 feeder 53-903 


24 feeder 53-1003 


25 factorty 103 
26 rect delta23 
27 rect wye23 


52 


0.000 100.000 0.0000 


0.000 
0.020 
0.060 
0.050 
0.050 
0.000 
0.050 
0.010 
0.020 
0.190 
0.020 


0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 


0.000 
0.000 


0.0000 


9999 
0004 


NDpweNdoOOIADMNAWNNY 


=OPlarnnhown 


a" 
tN 


26 15.35 543.78 0.00 0.00 
27 11.51 407.84 0.00 0.00 


9999 
0005 


007000001000020020 


0015 


0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 
0.00 0.000 


0.06600 30.00 
0.04100 0.00 
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0.000000 
0.000000 
0.000010 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
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